Bulk-edge correspondence in fractional Chern insulators 
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It has been recently realised that strong interactions in topological Bloch bands can stabilise novel states of 
matter. In this paper we study connections between these systems namely fractional Chern insulators, and the 
fractional quantum Hall states in the cylinder geometry using a generalised gauge-fixed Wannier-Qi basis. This 
new setup offers important advantages compared to the earlier exact diagonalisation studies on a torus. Most 
notably, it gives access to the properties of edge states and to the single-cut orbital entanglement spectrum, 
hence to the physics of bulk-edge correspondence. In addition, it is readily implemented in the state-of-the-art 
density matrix renormalisation group algorithm, which allows for numerical simulations of significantly larger 
systems. We demonstrate our general approach on examples of flat-band models on Ruby and Kagome lattices 
at bosonic filling fractions i/ = 1/2 and f = 1, where we observe the signatures of (non)-Abelian phases and 
establish the correspondence between the physics of edge states and the entanglement in the bulk. 
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Introduction. Fractional quantum Hall (FQH) states [1,2] 
provide some of the most unusual phases of matter. They sup- 
port excitations carrying only a fraction of the electron charge 
that obey anyonic statistics [3, 4]. One of the remarkable fea- 
tures of these states is that they realise a condensed matter 
example of the holographic principle [5], which in this con- 
text defines the relation between the physics of gapped bulk 
and of gapless edge states at the sample's boundary [2, 6]. 

Recently there has been a growing interest in the lattice 
models harbouring neaiTy-dispersionless (flat) bands charac- 
terised by non-zero Chern numbers, which can show states 
that are similar to ones observed in the fractional quantum 
Hall effect, so-called fractional Chern insulators (FCIs) [7- 
26]. These novel systems do not require external magnetic 
fields and can potentially be realised at room temperature due 
to shorter lattice length-scales compared to the typical mag- 
netic length in quantum Hall (QH) systems. While the con- 
nections between many-body correlations, quantum entangle- 
ment and the properties of edge states in the FQH case are 
relatively well understood [27-29], the coiTesponding physics 
in the FCIs has been relatively less studied. 

The standard numerical approach in the FCIs is the exact 
diagonalisation (ED) performed on small systems [10, 11]. 
Despite its considerable success in finding some of the ro- 
bust FCI states, there is a number of reasons to develop al- 
ternative numerical tools which would allow for a detailed 
understanding of larger systems. A potentially powerful 
method in this context is the density matrix renormalisation 
group (DMRG) [30]. This approach, which was initially de- 
signed for strongly-coiTelated one-dimensional systems, has 
been later successfully applied in the simulations of a variety 
of two-dimensional states of matter including geometrically 
frustrated magnets [31-33] and FQH systems [34—39]. 

From the recent developments in the theory of quantum en- 
tanglement, which brought new ideas of the area law [40], 



the entanglement spectrum [27], and the matrix product states 
[41^4], it has become clear that the cylinder geometry plays 
a very special role in the DMRG-based studies of strongly- 
coiTelated systems [41, 45]. Recently two pioneering papers 
suggested the application of DMRG in this geometry as the 
way to extract the information about topological properties of 
correlated states, including the FCIs [32, 46]. 

In this paper we consider a new setup which generalises 
the FCI description on a torus to the case of finite cylinders, 
by constructing the interaction matrix elements in the gauge- 
fixed version [47, 48] of the Wannier-Qi (WQ) basis [25]. Our 
approach has a number of advantages. First, it highlights the 
similarities with the standard FQH physics thus allowing for 
a direct comparison between the two [25, 47-49]. Second, it 
provides a computationally efficient setting for a momentum- 
space DMRG: for a lattice, whose linear dimension is N unit 
cells, the logarithm of the computational cost on a cylinder 
scales as N , compared with 2N for DMRG on a torus, and N"^ 
for ED [40]. Third, the orbital entanglement spectrum (OES) 
on cylinders probes the physics of a single edge, hence allow- 
ing for a cleaner and more straightforward identification of 
topological orders, compared to the torus setup which involves 
a non-trivial combination of two edges [50, 51]. Fourth, it nat- 
urally allows for the inclusion of an external potential which is 
ubiquitous in possible cold-atom realisations and, as we find, 
helps to stabilise some of the fragile FCIs including the coun- 
terpart of the non-Abelian Moore-Read state [52]. Finally, the 
presence of a physical boundary makes it possible to study the 
FCI edge states in details. 

Here we exploit all of these advantages by calculating the 
OES and the edge excitation spectrum for bosonic FCIs at fill- 
ing fractions v — 1/2 and v = 1. We find that the OES in 
these systems have the same low-lying counting structure as 
the edge excitation spectrum, thus providing a compelling ev- 
idence for the bulk-edge correspondence, similar to the one 



found in FQH states. In addition, we demonstrate that the 
Moore-Read FCI state at i/ = 1, which in a standard setup 
requires three-body interactions on a torus, is Hkely to survive 
for more realistic two-body interactions in the presence of a 
paraboHc potential, providing another possible way to realise 
non-Abelian phases which is different from the optical flux 
lattice setup of Refs. [18, 21]. We note that our calculations 
of the OES have been performed for systems which are much 
larger than the current limit of ED. 

Setup. We start with the lattice Hamiltonian on a finite torus 
whose periods are defined by two vectors Vi^2 with iVi.2 unit 
cells in the Vi.2 direction. The system consists of Nf, interact- 
ing bosons partially filling the lowest Bloch band of the Ruby 
[53] or the Kagome [7] lattice [we adopt the same symbols 
for the hopping parameters as used in Ref. [24] and fix their 
values to {U, U,tir, tii,U} = {1, 1.2, -1.2, 2.4, -1.46} for 
the Ruby lattice, and {ii, Ai, ^2, -^2} = {1,1,0,0} for the 
Kagome lattice]. We form a complete set of eigenstates in the 
lowest Chern band using gauge-fixed WQ orbitals localised 
in the Vi direction [25, 47]. The latter are counterparts of the 
lowest Landau level single-particle wave-functions in the Lan- 
dau gauge. This construction clarifies the connection between 
the lattice CIs and the QH systems. A generic translationally- 
invariant two-body lattice interaction (projected to the lowest 
band in a standard way [11]) has the following general form 
in the WQ basis 




Figure 1. (Color online) The ground-state orbital entanglement spec- 
trum for FCIs (red dots) and the corresponding OES for the Laughlin 
state (blue dashed lines) from DMRG for 18 bosons alu — 1/2. (a) 
The ground-state OES for 5 x 7 Ruby lattice is obtained for DMRG 
cutoff £ = 10~^" after 13 sweeps, (b) The ground-state OES for 
5x7 Kagome lattice (e = 10"^ after 10 sweeps). The OES of the 
Laughlin state is obtained with e = 10^^*^ after 20 sweeps. 
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where a' creates a boson in j's orbital. In order to ex- 
tend this construction to the cylinder geometry, we keep N2 
fixed while increasing A^i until we reach the convergence 
of torus matrix elements V,^'"'?' with j„ G [0, A^s — 1] for 
a given Ng <C NiN2- This simple procedure generates a 
TVj^' X iV2 lattice with Ng WQ orbitals on the finite-length 
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= Ns/N2- This system is described by 
which can be obtained from Eq. (1) by 
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iArj_j.oo ^s \ ■ Ths filling fraction v is defined as in terms 
of the number of particles, N^, and the number of "flux 
quanta", N^ = ^^^2, as i^ = Nb/{N, + S), where 5 is a 
size independent integer "shift". Compared with the Hamilto- 
nian of a FQH system on a finite cylinder the total momentum 
in our setup K — X]n=i i»i i^ conserved only mod iV2. 

Having constructed the Hamiltonian HH^ , we use DMRG 
and ED to study the FCIs numerically. Our DMRG imple- 
mentation is similar to the approach of Ref [37]. The orbital- 
cut entanglement spectrum of the ground state, which in the 
FQH case reflects the nature of edge excitations, can natu- 
rally be generated in finite-size DMRG sweeps. Using a stan- 
dard recipe developed for the FQH systems [54] we parti- 
tion WQ orbitals into two disjoint sets A and B, consisting 
of I A consecutive orbitals with "momentum" j running from 
to Z^ — 1 and the remaining Ng — I a orbitals with "mo- 



mentum" I A to Ng — 1. Generalized (mod N2) total momen- 
tum conservation requires that each OES level is labelled by 
Na and J a = iiT^ [mod iV2]. Below we will study Ruby 
and Kagome lattices of dimension N^^ x N2 on finite cylin- 
ders, whose QH counterpart is a cylinder with a circumference 
L = IbN2 •\/27r/ sin(7r/3), where Ib is the magnetic length. 

Entanglement spectrum. First, let us focus on the case of 
the filling fraction v — 1/2 in the presence of two-body on- 
site interactions X^i^-il"-* ~ 1)' where fii is the number of 
paiticles on i's lattice site. Taking into account that the stan- 
dard i^ = 1/2 bosonic Laughlin state on a cylinder appears at 
Ng = 2Nb — 1, we choose the same parameters for the FCI. 
Our ED calculation for H^^^ on small systems shows a unique 
ground state with excited states separated by the many -particle 
gap. The overlap |(^Gs|^Lau)| between the ground state and 
the Laughhn v = 1/2 state reaches 0.9996 and 0.9833 for 8 
bosons on 3 X 5 Ruby and Kagome lattices correspondingly. 
These strong overlaps signal that the ground states of HH^ are 
in the same class as the Laughlin state. 

To further identify the nature of topological orders in these 
systems, we use DMRG to calculate the OES for a cut in the 
WQ basis. By analogy with the FQH case the OES is ex- 
pected to reflect the properties of edge excitations. We choose 
I A = (^s + l)/2 and label each OES level by the particle 
number Na and the quasi-momentum J a of the A set. The 
ground-state OES is presented in Fig. 1 , where we have also 
included the results for the Laughlin state. Compared to the 
usual OES picture for FQH states, the one of the FCIs on a 
cylinder is folded, reflecting the generalised momentum con- 
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Figure 2. (Color online) The ground-state orbital entanglement spec- 
trum for FCIs in the Koo approximation (red dots) and the corre- 
sponding OES for the Laughlin state (blue lines) from DMRG for 25 
bosons at i/ = 1/2. (a) The ground-state OES in the A'oo approxi- 
mation for 7x7 Ruby lattice (e = 10~^^ after 9 sweeps), (b) The 
ground-state OES in the K^o approximation for 7x7 Kagome lat- 
tice (e — 10^^^ after 11 sweeps). The OES of the Laughlin state is 
obtained with s = 10^^" after 20 sweeps. 
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Figure 3. (Color online) (a,b) The ground-state orbital entangle- 
ment spectrum for FCIs in the Koo approximation for 16 bosons on 
3x5 Ruby lattice (red dots) with (a) three-body on-site interactions 
and (b) two-body on-site interactions in a confining potential with 
Vp — 0.01, compared with the OES of MR state (blue dashed lines) 
obtained from ED. In fact, even the finite-size reduced CFT-counting, 
indicated by numbers in parenthesis, match identically in the three- 
body case. (c,d) The ground-state OES in the K^o approximation for 
26 bosons on 5 x 5 Ruby lattice from DMRG for two-body on-site 
interactions with a confining potential Vp = 0.006 (e — 10~* after 
7 sweeps) in the (c) Na ~ 14, and (d) Na = 13 sectors. 



servation. This does not alter the counting of the low-lying 
spectrum, and one can clearly see the {1, 1, 2, 3, 5, 7} struc- 
ture (Fig. 1). The spectrum of the i/ — 1/2 FCI matches that 
of the corresponding Laughlin state with a very good preci- 
sion. 

Although the standard momentum K conservation does not 
hold exactly in the finite cylinder geometry, the ED results 
show that the FCI ground states still have a large weight 
(reaching 99.94% and 97.82% for 8 bosons on 3 x 5 Ruby 
and Kagome lattices) in the expected ii'-sector, i.e. K ~ 
Nb{Ni, — l) for a Laughlin state atu = 1/2. We use this obser- 
vation to justify the A'oo approximation, namely in the follow- 
ing we omit the terms that break the standard iiT-conservation 
in the Hamiltonian H^^^ . Then, we can obtain the ground 
states with a fixed K and label the OES by the standard 
{Na, Ka) pair (Fig. 2). One can see that the low-lying part 
of the K^ ground-state OES matches that of the ;/ — 1/2 
Laughlin state and has an entanglement gap (A^ sa 17 for the 
Ruby lattice and A^ « 13.5 for the Kagome lattice) which 
is larger than that for the i/ = 1/2 Coulomb ground state in 
the FQH on a sphere (A^ w 10 in this case [55]). These 
results, together with the folded ground-state OES provide a 
compelling evidence that the ground states of these systems 
at filling factor jy = 1/2 are in the same class as the FQH 
Laughlin state. 

Now let us consider the filling factor 1^ = 1 with Ng ~ Ni,— 
1, where the FCI counterpart of the non-Abelian Moore-Read 
(MR) state may potentially be stabilised. Compared with the 



Laughlin state, the MR-like phase in FCIs is more fragile, 
usually requiring three-body interactions on a torus [22-24]. 
Therefore, we first apply our cylinder setup to bosons with 
three-body on-site interactions J^i T^ii^T-i ~ l)(^i — 2), where 
fii is the number of particles on i's lattice site. The ED result 
for 10 bosons on 3 x 3 Ruby lattice shows a unique ground 
state with a large overlap |(^gs|*mr)| = 0.9997 with the 
exact MR state and the weight 99.96% in the expected K- 
sector, i.e. K = Nii{Ni,/2 — 1), which again justifies our K^ 
approximation. We also obtain remarkable ED results for a 
larger system: the overlap between the ground state in the 
Koo approximation and the exact MR state reaches 0.9998 
for 16 bosons and the low-lying part of the ground-state OES 
matches that of the exact MR state with very high precision 
and has a large entanglement gap A^ « 15 [Fig. 3(a)]. These 
results provide compelling evidence for the existence of FCIs 
with three-body interactions in the MR phase. 

Interestingly, we find that the MR FCIs might also be 
stabilised on cylinders by more realistic two-body on-site 
interactions in presence of a parabolic confining potential 
Vp T,^^o^ [2tt{j - ^^)/L]2a]aj in the bulk [52]. The over- 
lap between the ground state in the Koo approximation and the 
exact MR state reaches the value of 0.8829 for 16 bosons [Fig. 
3(b)], which is smaller than the three-body overlap, but is still 
quite high. The low-lying part of the OES has the {1,1,3} 
and {1,2} counting structures in two Na sectors, which is 
another signature of the MR phase, is also present in larger 
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Figure 4. (Color online) Edge excitation spectrum from ED in the 
Kac approximation aXv = 1/2 (two-body on-site interactions) and 
V = 1 (three-body onsite interactions), {a) Nb — 8 on 5 x 5 
Ruby lattice with vi = 0.01. (h) Nt = 8 on 5 x 5 Ruby lat- 
tice with v'p = 0.0006. (c) Nb = 8 on 5 x 5 Kagome lattice 
with v'p = 0.001. We observe the U(l) counting of low-energy 
excitations {1, 1, 2, 3, 5, 7, ...} at AK = {0, 1, 2, 3, 4, 5, ...} in all 
three cases, (d) A^;, — 10 on 5 x 3 Ruby lattice with v'p = 
0.0002. We observe the U(l)xlsing counting of low-energy exci- 
tations {1, 1, 3, 5, 10, 16, ...} at AK = {0, 1, 2, 3, 4, 5, ...}. 



systems which we study using DMRG, see Fig. 3(c,d). 

Edge excitation spectrum. The open boundaries on finite 
cylinders provide a natural setting for the studies of edge ex- 
citations which appear in the vicinity of real physical edges. 
In our numerical approach we keep N2 fixed and increase 
N'^ adding extra WQ orbitals, then open the edge on one 
side to allow occupation of these states while keeping the 
edge on the other side closed. In order to observe a stable 
edge excitation spectrum, we consider various confining po- 
tentials that extend from the bulk to the extra WQ orbitals. 
For not too strong potentials, a branch of low-energy exci- 
tations separated from higher levels appears in the spectrum 
for both filling fractions v = 1/2 and j/ = 1 as shown in 
Fig. 4. For a linear confinement, vi X]j=o i'^''^J/L)ajaj, the 
spectrum accurately matches the prediction of Luttinger liq- 
uid theory [6]: the dispersion is linear and the edge states in 
each AK sector have nearly-degenerate energies [Fig. 4(a)]. 
This degeneracy can be split by a parabolic confinement, 
^p^i=o (27rj7L)^aJaj, which makes the excitation spec- 
trum similar to the OES [Fig. 4(b-d)]. The number of edge 
states in each AK sector does not depend on the form of the 
confinement and matches exactly with the conformal field the- 
ory prediction until the finite-size effects intervene at higher 
energies. Similar results have recently been obtained for a re- 
lated problem of FQH states on a lattice in uniform magnetic 
field [56]. 

Discussion. In summary, we studied bosonic fractional 
Chern insulators in the finite cylinder geometry using a combi- 
nation of the exact diagonalisation and the momentum-space 



DMRG. The ground-state OES at 1/ = 1/2 has a strong over- 
lap with the OES of the corresponding Laughlin state. The 
ground-state OES at i^ = 1 shows that the FCI analogue 
of the FQH Moore-Read state is likely to survive even with 
two-body on-site interactions. The counting structure in the 
ground-state orbital entanglement, and the edge excitation 
spectrum, provide strong evidence for the bulk-edge coiTe- 
spondence in FCIs. Our setup is likely to bring new insights 
into intriguing and less understood FCI states which have no 
direct QH counterparts, most notably the states which can ex- 
ist in flat bands with higher Chern numbers [57-59]. 
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